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We suggest how we can use the mass profile of galaxy clusters beyond their virial radius to 
measure their mass accretion rate, a key prediction of structure formation models. The mass 
profile can be estimated by applying the caustic technique to dense redshift surveys of clusters 
and their outskirts, where dynamical equilibrium does not necessarily hold. An additional probe 
of the mass growth of clusters is their mass fraction in substructures. We show that the caustic 
technique, that identifies cluster substructures as a by-product, returns catalogs of substructures 
with mass larger than a few \Q l 'p l M@ that are between 60% and 80% complete, depending on 
the density of the redshift survey. 
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Figure 1: Current measures of f{z)o g at different redshifts. The thick solid (dashed) curve shows the 
General Relativity prediction in a ACDM model with WMAP9 (Planck) parameters, while the dotted, dot- 
dashed, and dot-dot-dashed curves show the DGP, coupled dark energy, and f(R) models, respectively. From 
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1. The mass accretion rate of galaxy clusters 

In the next decade or two, ongoing and upcoming wide-held imaging and spectroscopic red- 
shift surveys (e.g., DES, eBOSS, DESI, PFS, LSST, Euclid, WFIRST) aim to measure the growth 
rate of cosmic structure on linear and mildly non-linear scales, up to wave numbers k ~ 0.2 h M pc 1 , 
in the redshift range 0 < z <2. A quantity that has been commonly measured is f(z)o%, where 
f(z) = dlnD/dlna, a = 1/(1 + z) is the scale factor, D(a) the linear growth factor, and the 
normalization of the power spectrum of the density perturbations. The claimed accuracy of the 
growth rate measured, for example, with Euclid is 1% to 2.5% in the redshift range 0 < z < 2 [[I]]; 
however, current measures of f(z)Os based on redshift space distortions up to redshift z = 0.8 are 
affected by uncertainties between 10% and 50% |[|]] (Figure |IJ). Combining galaxy-galaxy lens- 
ing with redshift-space distortions is a promising improvement that, however, is not yet accurate 
enough to discriminate among different cosmological models and modified gravity theories [[|| . 

On very non-linear scales, above the wave numbers k ~ 0.2 li Mpc 1 . the measurement of the 
growth rate of individual dark matter halos has only been attempted on the scales of galaxies: by 
combining the number of observed pairs of close or disturbed galaxies with the merger probability 
and time scale derived from /V-body simulations, we can infer the galaxy-galaxy merger rate [ 5|. |], 
0. However, current results are inconclusive |^j|: in fact, the merger rate of dark matter haloes and 
the merger rate of galaxies are related by dissipative processes that are difficult to model [|j] and 
these two rates do not necessarily coincide 00 Ell- 

Dissipative processes are less relevant during the mass accretion of galaxy clusters, whose 
rate could be simply estimated based on the measurement of the amount of mass in the cluster 
outskirts. However, this advantage over the galaxy-galaxy merger rate has never been capitalized 
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because, in the large and less dense cluster outskirts, (1) cluster galaxies are difficult to distinguish 
from foreground and background galaxies, and (2) other probes, e.g. X-ray emission, are below the 
sensitivity of current instruments. Moreover, the cluster outskirts are not in dynamical equilibrium 
and the usual mass estimation methods based on virial equilibrium are inappropriate. 

This observational deficiency clashes with the numerous and detailed studies of the mass 
growth of galaxy clusters in V-body simulations based on the identification of their merger trees 
[|, 0, 0, 0, 0], Various laws for the average mass accretion history of dark matter halos 
have been proposed: (1) M(z) = Moexp(— az) where Mo/2 is the halo mass at the formation red- 
shift of the halo Zf = \n2/a 0; (2) log[M(z)/M 0 ] = -0.301 [log( 1 + z)/log(l +Zf)] n l[f2|l; (3) 
M(z) = M 0 (l +z) b exp(-rz) [0 0], The discrepancies between the various relations are due to 
different mass and time resolutions and different halo statistics of the V-body simulations used to 
infer the relations. 

The Caustic group 1 in Torino has started a project aiming at estimating the mass accretion 
rate of galaxy clusters by measuring the mass of a spherical shell surrounding the cluster and its 
infall time. This approach is rather crude when compared with the stochastic aggregation of dark 
matter halos in the hierarchical clustering formation models. Nevertheless, our preliminary results 
are promising: they suggest that measuring the mass accretion rate of galaxy clusters is actually 
feasible and can potentially provide a new observational test of the cosmological and structure 
formation models. 

2. Cluster mass profiles in the outer regions: the caustic method 

If we can measure the mass profile of a cluster beyond its virial radius r?oo, we can estimate 
its instantaneous mass accretion rate as 

M = M[< r20o(1 + 5j)] ~ M2Q0 (1 +z) 3 / 2 (2.1) 

where M 200 = M(< /too). Equation ( |2. l| ) assumes that a shell of proper radii r 2 oo and r 2 oo(l + 5?) 
takes a cosmic time At to fall onto the cluster with constant acceleration — GM 200 /KI + 5,/2)r2oo] 2 
and null initial velocity. We keep At = 0.1 Gyr fixed with redshift and accordingly vary S s = 
100H 2 (z)A t 2 , where H(z) is the Hubble parameter. The additional factor (1 + z) 3/2 corrects for 
the transformation from the cosmic time At to the infall time derived with the proper radii r 2 oo and 

>-200(1 + 5,) 0. 

Figure [2] shows how this simple recipe compares with the CoDECS V-body simulation [ |~iT> j| 
of clusters in a ACDM model. The thick line shows the mean M computed with equation ( fP| ) 
for a sample of clusters with M 200 = 1 0 14 h 1 Mr., at z = 0, whereas the thin line is the mean mass 
accretion rate obtained from the merger trees of the dark matter halos. The mean accretion rate from 
the merger trees lie in the region of one standard deviation of our M (shaded area). In addition, the 
mean rate from the merger trees is recovered by the mean of our prescription within 20% in the 
redshift range z= [0,1]. 

The results of Figure |] show that estimating the mass accretion rate of clusters is indeed fea¬ 
sible if we can measure the cluster mass profile beyond /too- This measurement can be performed 

*www.dfg.unito.it/ricerca/caustic 
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Figure 2: Comparison between our estimate of the mass accretion rate of a sample of M 200 = 10 14 /i 1 AC. 
clusters and /V-body simulations. The shaded area shows the 1 a range around the mean mass accretion rate 
estimated with equation (2.1) (thick solid line). The thin solid line is the mean mass accretion rate derived 
from the merger trees. Adapted from JTsIl . 


with the caustic technique. The caustic technique [|20|, |2J_[ [22|] uses the celestial coordinates and 
redshifts of galaxies to estimate the gravitational potential and mass profiles of a cluster from the 
central region to radii much larger than the virial radius. The caustic technique (i) does not rely on 
the dynamical equilibrium of the cluster; (ii) measures the three-dimensional distribution of mass, 
based on the assumption of spherical symmetry; (iii) estimates a mass that is unaffected by the 
presence of substructures within the cluster | fH] |22| | and by structures along the line of sight [23]; 
(iv) can be applied to clusters at any redshift and it is only limited by the telescope time required to 
measure a sufficiently large number of galaxy redshifts. A robust estimate of the cluster mass out 
to 3r2oo requires ~ 200 galaxy redshifts [[l2||. Figure || shows how, with this number of redshifts 
in the field of view of a cluster, the caustic technique recovers the mass profile up to 4/-2oo with 
no bias (solid squares) and with ala relative uncertainty of 20% (error bars). The required num¬ 
ber of galaxies for such an accurate estimate was demandingly large in the late nineties, when the 
technique was designed, but, thanks to the development of multi-fiber spectroscopy, it is a feasible 
target nowadays. 

The measurement of a cluster mass accretion rate is tightly linked to the estimate of the total 
mass of the cluster within its turnaround radius, the so-called ultimate mass M u [^5|]: with the caus¬ 
tic technique, by combining the 50 CIRS clusters [25] with the 58 FleCS clusters [26], we found 
M u /M 200 = 1.99 ± 0.11, a measure accurate to 5% [|26|]. Our measure agrees with the ACDM pre¬ 
diction, where M u /M 2 oo has a log-normal distribution with a peak at mass ratio 2.2 and dispersion 
0.38 [27]. The accuracy of this unique estimate of the ultimate mass M u suggests that we can aim 
to measure the mass accretion rate with the caustic technique to a 20% accuracy or better. 

In principle, we could also estimate the mass in the outskirts of clusters with weak gravitational 
lensing analyses. However, the major source of uncertainty with this approach is the mass projected 
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Figure 3: Mass accuracy with the caustic technique: the figure shows the average ratio between the caustic 
estimate based on 200 galaxy redshifts and the true mass of simulated clusters. Error bars show the lcr 
uncertainties that mostly originate from projection effects. The overestimate at r < rioo is due to a known 
systematic error [[22]]: at any rate, it is irrelevant in the estimate of the accretion rate with equation (2.1), that 
involves a mass difference. The three error bars at r ~ r 2 oo show the uncertainty on M 200 estimated with 
weak lensing for a 10 I5 /i _ 1 M@ cluster at z = 0.3 (20%), z = 0.1 (50%) and for a 5 x 10 14 /z _1 M 0 cluster at 
z = 0.1 (80%) 
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Figure 4: Comparison of the weak lensing mass profile for MS0906 with 1 a errors (heavy and light dashed 
curves respectively, [[28]]) with the mass profiles derived with the caustic technique for MS0906 (open di¬ 
amonds), A750 (open squares), and the effective sum of the caustic mass profiles for the two superposed 
clusters taking the 0.6/z 1 Mpc offset between the centers into account (solid squares). From [[23||. 
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along the line of sight. This uncertainty increases with decreasing cluster mass and increases with 
increasing radius. More importantly, the uncertainty increases at both lower and higher redshifts. 
Figure |3] shows that the expected accuracy of the caustic technique is a factor of two better than the 
accuracy expected with weak lensing for clusters at z = 0.3, and a factor of four better for clusters 


at z = 0.1 | [24[ . [29| ]. Finally, for the accretion rate estimate, the caustic method can be applied to 
clusters at any redshift, whereas weak lensing is basically limited to clusters within the redshift 
range 0.1-0.5 


,29]. 


Figure [| shows a remarkable example of how a superposition along the line of sight can affect 
the weak lensing mass: the two clusters A750 and MS0906 are almost aligned along the line of 
sight and the weak lensing mass of MS0906 is approximately the sum of the masses of the two 
clusters; a sufficiently dense redshift survey shows that the two clusters are distinct and the caustic 


technique can easily estimate their individual masses [23]. 


3. Identification of cluster substructures with the caustic technique 

An additional diagnostic of the mass growth of clusters is the presence of substructures in the 
central region and the outskirts of the cluster. The first most serious difficulty of this approach is 
the identification of substructures. Suggested methods are based on the distribution of galaxies in 


redshift space [[30||, X-ray emission [|3T], [32||, and gravitational lensing [ [33| | (see [ |34| | for additional 
references). Some investigations have attempted to use this piece of information to constrain the 
dynamical state of clusters or their formation time (see e.g., [| 


, 36 


)• 


With optical observations, the best approach clearly is to combine galaxy positions and red- 
shifts. The first step of the caustic method is to arrange the galaxies in a binary tree based on 
a galaxy pairwise projected binding energy, similarly to the method proposed by Serna and Ger- 
bal [|39[]. However, the caustic method goes further and identifies two thresholds that separate the 
branches of the binary tree and provide a list of groups and a list of substructures from the cluster 
center out to its outskirts. 

The substructures identified in redshift space can correspond to real substructures or can be 
due to chance alignment; in addition, some of the real substructures of the cluster may not be 
identified. The number of false detections and the completeness of the substructure sample depend 
on the cluster mass, on the substructure mass and, crucially, on the density of the redshift survey. 
When applied to mock redshift surverys of clusters extracted from X-body simulations, with typical 
properties of surverys like CIRS [25] and HeCS [26], namely ~ 200 redshifts within 3r2oo and 
cluster mass M 200 r^j 10 14 h~ l M Q , r^j 50% of the substructures identified with the caustic technique 
are false detections, while the completeness is between 60% and 80% for substructures masses 
larger than a few 10 13 /; 'AT. [|34]| . 

We thus conclude that the caustic technique appears to be a very promising method to identify 
substructures in galaxy clusters. In future work, we will investigate how this technique compares 
with methods based on X-ray and gravitational lensing studies. 


I am grateful to my collaborators Ana Laura Serna, Cristiano De Boni, and Heng Yu for making 
this old idea come true and for allowing me to show results from our common projects. I thank 
Christian Marinoni and the other organizers of Frontiers of Fundamental Physics 14 for giving me 
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